Structure of many-body continuum states in the proximity of exceptional points 
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We demonstrate existence of exceptional points in many-body scattering continuum of atomic 
nucleus and discuss their salient effects on the example of one-nucleon spectroscopic factors. 
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INTRODUCTION 

Degeneracies of eigenvalues of matrices dependent on 
parameters is a problem of great importance in Physics. 
The interest in this problem goes back to Hamilton 
who pointed out a conical refraction due to coincident 
eigenvalues [2]. Later, the significance of degeneracies 
of matrix spectra has been revealed in quantum physics 
B S]' crystal optics Q and chemistry 0]. Eigenvalue 
degeneracies are also essential for an understanding of 
the spectral fluctuations 0] and the onset of quantum 
chaos Q • Among those degeneracies, one finds the excep- 
tional points (EPs) [1, [^, where two Riemann sheets 
of the eigenvalues are entangled by the square-root type 
of singularity. The EPs are generic singularities and, 
hence, can be found in numerous situations, e.g. the 
resonance problems 11, l2| or the quantum phase tran- 
sitions 



3,111, El- 
Much effort has been devoted to studies of degenera- 
cies associated with avoided crossings in quantal spectra, 
focusing mainly on the topological structure of Hilbert 



space and the geometric phase [11|, |16|, [17[ which a 
quantum system acquires when transported adiabatically 
around the singularity in parameter space 3, 3 ■ These 
novel effects have attracted a considerable experimental 
interest as well fl^ . 

In many-body systems, EPs have been studied in 
schematic models, such as the Lipkin model the 
interacting boson model 1J| and, recently, the 3-level 



pairing model [21|, |24| which belongs to a general class 
of the Richardson-Gaudin models [23l [23 |. In the lat- 
ter case, it has been found that eigenvalue crossings in 
the complex-extended parameter space of the Hamilto- 
nian and the reduction in the number of EPs uniquely 
define the quantum integrable system [2l|. If true, this 



conjecture opens a possibility for studies of fingerprints 
of the chaotic dynamics in the quantum regime for small 
dimensional Hilbert spaces. 

In the case of closed quantum system (CQS) described 
by a hermitian Hamiltonian, mixing of two eigenvalues 
leads generically to the level repulsion and the avoided 



level crossing. In the open quantum system (OQS), 
eigenvalues corresponding to unbound levels may exhibit 
not only crossing of either energies or widths but also 
true degeneracies due to a joint crossing of both ener- 
gies and widths. The latter case corresponds to a de- 
generate double pole of the scattering matrix (S-matrix) 
Wl . 27, 28, 2^. Hence, OQS is a natural system to look 
for singularities of the EP type. 

First attempt to search for EPs in the realistic many- 
body model of atomic nucleus has been made recently 
using the real-energy continuum shell model [sOl , the so- 
called Shell Model Embedded in the Continuum (SMEC) 



31, 32 



33|, which is defined in projected Hilbert spaces 
with at most two particles in the scattering continuum. 
The total function space of the A-particle system in 
SMEC consists of two sets: the set of square-integrable 
functions {'I'f }, used in the standard nuclear Shell Model 
(SM), and the set of continuum states (the environment 
of the CQS) {Cb"^-*} embedding the SM. In one-particle 
continuum problems, channels c are determined by the 
motion of an unbound particle in a state Ij relative to 
the A — 1 nucleus in a certain SM state ^f~^- An in- 
dex denotes an outgoing asymptotics of a continuum 
state . Eigenstates of the OQS are found by solving 
the eigenvalue problem for an effective Hamiltonian: 



(1) 



in the function space {^f} of discrete states. Indices 
Q and P in ([1]) denote the projected Hilbert subspaces 
with and 1 nucleon in the scattering continuum, respec- 
tively. G^'' (E) is the Green function of a single nucleon 
in the P subspace. Whereas the SM Hamiltonian {Hqq) 
is hermitian, Hqq is an energy dependent, non- hermitian 
operator above the particle emission threshold and a her- 
mitian operator below the emission threshold. For each 
total angular momentum J and parity tt separately, the 
Q — P coupling term in ([1]) gives rise to an external mix- 
ing of different SM states {^ki'^^)} coupling to 
common decay channels, both closed and opened. The 
amount of this external mixing depends on the position 



2 



of considered SM states with respect to the decay channel 
thresholds [3J]. More details about the SMEC approach 



and its recent apphcations can be found in Refs. [25ll3 
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FIG. 1: Energy and width of lowest two eigenstates Oj , Oj of 
an effective Hamiltonian ([T]) in ^®Ne. These eigenstates form 
an exceptional point at i? = 26.523 MeV. The continuum- 
coupling strength corresponding to this degeneracy is Vb = 
-1182.32 MeV-fm^. Energies (in MeV) on the vertical Re_E,- 
axis are plotted with respect to an arbitrarily chosen origin. 
Ifence, only relative values of the real part of energies have a 
meaning in this representation. On the contrary, the width 
Ti of each eigenstate depends only on a relative distance from 
the elastic threshold {E = 0) and is not influenced by the 
convention chosen to plot Rei5i. 



In this work, we study the one-proton continuum of 
^^Ne to look for a possible EP with = 0+ in a phys- 
ically relevant region of excitation energies and param- 
eters of the SMEC. For the SM Hamiltonian {Hqq in 
(P)) we take the ZBM Hamiltonian [35]. The resid- 
ual coupling Hqp between Q-subspace and the embed- 
ding continuum is given by: V12 — Vb(5(ri — r2). The 
true degeneracy of two 0^ resonances is searched for 
by varying the strength parameter Vq in Hqp and the 
total energy E of ^^Ne. i? = is chosen at the low- 
est one-proton emission threshold. In present studies, 
we consider the coupling of 0+ SM states to three low- 

r 

est decay channels in ^^Ne: '^^'F{1/2'1) ® P£;(lsi/2) 



(the elastic channel), ^^F(5/2J') ® p_E(0d5/2) 



and 



15 



F(l/2r) ® P£(0pi/2 



eff 



The energy dependence of 

16^ 



QQ is due to the coupling of discrete states of °Ne 

15t 



to the decay channels 



This 
eff 



coupling depends explicitly on E. For E > Q, Hqq is 
complex-symmetric and its eigenvalues are complex. 
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FIG. 2: Real and imaginary parts of the one-proton spec- 
troscopic factor (^^Ne(0+)|[^'^F(l/2+) ® p(lsi/2)]^> to the 
ground state {.F = 1/2+) of ^^F in diflterent states 



,4) of '"Ne (Fo 



-1100 MeV-fm'' 



O2 eigen- 



Fig. [T] shows energies and widths of 0^, 
values as a function of the total energy of ^^Ne for a 
fixed value of the strength Vq = -1182.32 MeV-fm^ 
of the residual couphng Hqp. Oj^, 0^ eigenstates and 
eigenvectors of H^{E) coalesce at E = 26.523 MeV, 
forming an EP. Encircling this singularity in a parame- 
ter space of Hqq, reveals a topological phase in 0^, 
eigenvectors. One can see a strong variation of 0^ and 
O2 eigenvalues in the vicinity of thresholds for channels 

0+ r 1O+ 
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F(l/2+)®p(lsi/2) 



and 
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F(5/2+)®p(0d5/2) 
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FIG. 3: Real and imaginary parts of the one-proton spectro- 
scopic factor (^®Ne(0+)|[^'5F(5/2f ) (g)p(0d5/2)]"^> to the first 
excited state (J" = 5/2^'') of ^^F in different 0+ states 
(n = I,-- - ,4) of i^Ne {Vo = -1100 MeV-fm^). 
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FIG. 4: Real and imaginary parts of the one-proton spectro- 
scopic factor (^'^Ne(0+)|[^^F(l/2j") ®p(0pi/2)l^) to the sec- 
ond excited state (J'' = 1/2^) of ^'^F in different 0+ 
(n = 1, • ■ ■ , 4) of I'^Ne (Vo = -1100 MeV-fm^). 



states 



respectively. The Coulomb interaction shifts the near- 
threshold minima of RcE'i {E) slightly above the thresh- 



old in channels 



[E ^ 0) for 0+ 



F(l/2+)®p(lsi/2) 
15f(5/2+) ® p(0(i5/2) {E = 0.671 MeV) for 



and 

One should also notice that the energy E of the lowest 
EP with = 0+ does not necessarily correspond to the 
resonance energy which in SMEC is determined by solv- 
ing a fixed-point equation for eigenvalues Ei{E) of the 
effective Hamiltonian Hqq{E) [32|. Real part of ener- 
gies in Fig. [T] are plotted with respect to an arbitrarily 
chosen zero on the axis ReE'i and only relative differences 
are meaningful. In the case considered in Fig. [U if the 
ground state energy with respect to the elastic threshold 
{E = 0) is fixed at £^0+ = -0.5 MeV or -^0.05 MeV, 

then the energy of the 0^ state obtained from the fixed- 
point equation becomes i?Q+ — 4.918 MeV or 6.251 MeV, 
respectively. 

One should stress an essential role of non-hermiticity 



of E[qq in this problem. The hermitian Hamiltonian can 
be always diagonalized and eigenvalues, even if degener- 
ate, always correspond to distinct eigenvectors. On the 
contrary, in the non-hermitian Hamiltonian one may find 
non-trivial Jordan blocks of eigenfunctions and singular 
points in a parameter space of the Hamiltonian where 
both eigenvalues and eigenvectors coalesce. These EPs 
have drastic effects on systems behavior, as we shall ex- 
amine below on the example of spectroscopic factors. 

Fig.[5]shows real and imaginary parts of the one-proton 
spectroscopic factor (i6Ne(0+)|[i5F(l/2+) ® p(lsi/2)]^) 
for the four lowest 0+ states in ^^Ne. This spec- 
troscopic factor comes from a coupling of 0+ states 
(n = !,•■•, 4) of I'^Ne to the ground state 1/2^^ of 
^^F. The strength of the residual (5-force coupling be- 
tween Q-subspace and the embedding continuum (Vb = 
— 1100 McV-fm'^) is chosen to reproduce an experimen- 
tal width of 5/2j^ state in ^^F. Analogously, Fig. [3] 
and HI exhibit the (^^NeiO+)\[^^F{5/2+) (g) p{0d5/2)f ) 



and 



/16t 



Ne(0+)|rF(l/2J-)(8)p(0pi/2)] > spectroscopic 
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factors, respectively. These spectroscopic factors result 
from a coupling to the first (J'^ = 5/2^) and the second 
{J'^ = 1/2 j") excited states in the daughter nucleus -^^F. 
Imaginary part of the spectroscopic factor, which appears 
for energies above the first one-proton emission threshold 
(i? > 0), yields the uncertainty of a real part and is re- 
lated to the decaying nature of considered 0+ states. For 
deeply bound states, an external mixing of SM states due 
to the coupling to common decay channels can be safely 
neglected and, hence, the spectroscopic factor calculated 
in SMEC becomes close to the SM value. This tendency 
can be seen in Figs. [2]|4]for £' < 0. 

Figs. [2ll4l exhibit the external mixing of = 0+ SM 
states by showing the i?-dependence of spectroscopic fac- 
tors for corresponding 0+ SMEC eigenstates. This mix- 
ing is strongly dependent on the structure of individual 
0+ SM states, the position of different one-proton emis- 
sion thresholds, and the total energy E of ^^Ne. One 
can see in Figs. [21 [3] that the one-proton spectroscopic 
factors are large for 0^ and 0^. The eigenstate Oj^ has 

only a sizeable {^^^e{Q+)\[^^Y{l/2^) ®Y>{'dp-^/2)V) spec- 
troscopic factor. 

One can notice two distinct phenomena in the in- 
dependence of spectroscopic factors: the threshold ef- 
fect and the EP effect. The particle-emission threshold 
is a branch point in OQSs [3a]. In Figs. [2]|4l one may 
notice a strong effect of the coupling to all three decay 



channels 



i5F(l/2+). 



15 



F(5/2+)®p(0d5/2) 



P(lSl/2^ 

0+ 

, and 



(the elastic channel), 

i5F(l/2r)®p(0pi/2)]° 
Remaining O"*" states 



for a pair of eigenstates 0^, OJ 
are effectively decoupled showing a weak i^-dependence 
of spectroscopic factors in the vicinity of thresholds. The 
threshold effect is well localized in E and leads to mutu- 
ally compensating variations in spectroscopic factors for 
strongly coupled 0^ and 0^ eigenstates of These 
compensating variations of single-particle occupancies is 
a signature of an entanglement of the two externally cou- 
pled states in the scattering continuum. In other words, 
one cannot define any matrix element or expectation 
value of an operator which does not commute with the 
Hamiltonian on an individual state of the entangled cou- 
ple of eigenstates. 

The effect of an EP at = 26.523 MeV leads 



spec- 



rNe(O+)|[i5F{5/2+)0p(Od5/2r), 
C«Ne(0+)|[i5F(l/2r)®p(0pi/2)]") 



to much stronger variations of all three 
troscpic factors: (if^Nc(0+)|PF(l/2+) » p(lsi/2)]"), 

and 

for 0+ and 0^ 
states. The same compensation mechanism as seen in 
the vicinity of the thresholds is found for spectroscopic 
factors in 0^ and states. However, since the strength 
of the residual couphng (Vb = —1100 MeV-fm^) does not 
correspond to an exact value (Vb = —1182.32 MeV-fm'^) 
for which the true degeneracy of eigenvalues has been 



found, the variations seen are caused by the proximity 
of the EP and not by the EP itself. As a result, the EP 
effect is weaker in magnitude and less concentrated in 
energy E than expected at an exact position of the EP. 
It is interesting to notice that the intricate entanglement 
of two continuum states caused by the EP is seen 
not only at the singularity point in parameter space 
of the Hamiltonian {E = 26.523 MeV, = -1182.32 
MeV-fm'^) but also relatively far away from it. In that 
sense, it is a robust phenomenon relatively weakly de- 
pendent on a fine-tuning of the Hamiltonian parameters 
and apt to be observed in nuclear physics experiments. 

We acknowledge fruitful discussions with J. Dukelsky 
and W. Nazarewicz. 
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